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Abstract Wc apply a recently-developed coarse-grained 
model of DNA, designed to capture the basic physics 
of nanotechnological DNA systems, to the study of a 
'burnt-bridges' DNA motor consisting of a single-stranded 
cargo that steps processively along a track of single- 
stranded stators. We demonstrate that the model is 
able to simulate such a system, and investigate the sen- 
sitivity of the stepping process to the spatial separa- 
tion of stators, finding that an increased distance can 
suppress successful steps due to the build up of un- 
favourable tension. The mechanism of suppression sug- 
gests that varying the distance between stators could 
be used as a method for improving signal-to-noise ra- 
tios for motors that are required to make a decision at 
a junction of stators. 
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1 Introduction 

DNA has proved itself to be a versatile material for en- 
gineering on the nanoscale. A single DNA strand con- 
sists of nucleotides covalently connected by a sugar- 
phosphate backbone. Nucleotides also possess one of 
four types of base group, adenine (A), guanine (G), cy- 
tosine (C) or thymine (T). The planar bases interact 
with other bases by forming parallel face-to-facc stacks 
and hydrogen bonds along their edges, resulting in the 
iconic double-stranded DNA helix (dsDNA) [1]. Due to 
their chemical structure, hydrogen-bonded C-G and A- 
T base pairs are the most favourable, and helices made 
from complementary sequences of bases form the most 
stable duplexes. 

The selectivity of DNA base-pairing can be used to 
engineer artificial structures and devices. As was origi- 
nally proposed by Seeman, sequences of a set of single 
strands can be designed so that a certain configura- 
tion is the global free-energy minimum of the system 
[2]. An enormous range of nanostructures have been 
realised simply by cooling solutions of single-stranded 
DNA (ssDNA) based on this principle. Finite-size struc- 
tures have been designed that assemble from a small 
number of short oligonucleotides [3| , a single long strand 
and many short staple strands [4l[5] (a technique known 
as DNA 'origami') and a large number of short stands 
[51[7]- Additionally, large one-dimensional ribbons [5], 
2-dimensional arrays [^[TU] and 3-dimensional crystals 
[TT] have been realised. 

DNA nanotechnology is not confined to static struc- 
tures - dynamic systems have also been realised. Du- 
plex formation and toehold-mediated strand displace- 
ment |12| (a process in which a strand is removed from a 
duplex by a competing strand that can form more base 
pairs with the complement) allow a DNA system to re- 
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spond to its environment. In particular, these processes 
can couple chemical change to mechanical operations 
and have the potential to process signals. DNA nan- 
otweezers [13], a switch that can be cycled through its 
closed and open states by the sequential addition of two 
types of strand, demonstrated the principle. "Clocked" 
addition of strands P^fBlfTCllTTl lTS l fT^ or permutation 
of external conditions pOll^ have since been used in 
the design of a number of active systems, including 
some in which the mechanical change has been har- 
nessed to induce unidirectional motion along a track 
[T8l[T9l[2T] . Recently, autonomous devices and walkers 
that function without external forcing (by catalysing 
the equilibration of an out-of-equilibrium system) have 
also been created Pl[^[^[^[^[?7l[^[^|5ni|5T] . 

The high degree of parallelisation and the ability to 
interface directly with biological and molecular systems 
make DNA-based computation promising. In 1994, Adle- 
man showed that DNA strands could be used to encode 
a Hamiltonian path problem, which was then solved 
upon mixing of the strands [22 • Since then, much work 
has gone into developing DNA-based logic circuits [33] , 
with a DNA neural network that can recognise sim- 
ple patterns having recently been developed [31] . DNA 
logic has also been combined with walking devices to 
produce systems that can select from distinct pathways 
at a junction depending on solution conditions, or prop- 
erties of the walker itself j3 11155] . 

Although DNA nanotechnology shows great promise 
as a field, it is far from being fully mature. In particular, 
optimisation of nanostructure assembly and nanodevice 
operation will be vital if they are to prove generally use- 
ful. DNA nanostructures and nanodevices are typically 
designed using well-established thermodynamic models 
of DNA duplex stability, such as the unified nearest- 
neighbour model of SantaLucia }36| . However, nanode- 
vices and nanostructures can involve non-trivial multi- 
stranded complexes with pseudoknots |37| or complex 
internal loops whose stabilities have not yet been incor- 
porated into thermodynamic models. Moreover, non- 
equilibrium processes can be important in these sys- 
tems. Furthermore, the three-dimensional structure of 
a DNA complex may result in tension or compression 
forces [55] that cannot be described without an explicit 
three-dimensional representation of the system. 

Computer simulations provide controllable access 
to time and spatial resolutions that are not accessi- 
ble in experiments. Simulations of DNA nanotechnolog- 
ical systems therefore have the potential to offer insight 
into aspects of their operation and design, provided the 
computational model accurately describes the relevant 
physics. At the greatest level of detail, DNA molecules 
can be studied using quantum-mechanical calculations. 



However, typical of what is achievable with the most 
advanced methods is for example the calculation of in- 
teractions for two base pairs in vacuum |39| . Classical 
molecular dynamics simulation packages such as AM- 
BER [40] retain an all-atom representation of the DNA 
molecule and use empirical force fields between atoms. 
However, they too are computationally very demand- 
ing, and simulating rare events, such as breaking of 
base pairs in a duplex, remains at the limits of what 
is currently possible jH]. 

Processes relevant to DNA nanotechnology typically 
involve reactions between multiple strands whose lengths 
vary from about ten to several thousands of bases. In 
order to address sizes and time scales at which the re- 
actions relevant to DNA nanotechnology happen, one 
needs to turn to coarse-grained models of DNA [^H51 
|iIl|i51|i51|T7l|i51|i^[5gi[n][5^[551[5il[55][551[57] . Such mod- 
els replace multiple atoms by a reduced set of degrees of 
freedom and parametrised effective interactions. Each 
model necessarily makes a compromise between the level 
of detail of its representation, physical accuracy and its 
efficiency, which sets its specific domain of applicability. 

In this paper we use the coarse-grained DNA model 
of Ouldridge et al. [55115^150] to study the operation of a 
'burnt-bridges' DNA motor created by the Turberfield 
group [M1I5II1I55] . The model, which provides a physi- 
cally reasonable representation of ssDNA, dsDNA and 
the transition between the two, allows for simulations 
involving reactions of DNA systems of order tens to 
hundreds of bases in total on a single CPU and several 
thousands on a GPU. Further, it has been previously 
successfully applied to the study of DNA nanotechno- 
logical systems including DNA tweezers [HI], kissing 
hairpins [SI] and a two-footed DNA walker [SHj . 

We first give a brief overview of the design of the 
burnt-bridges motor and the DNA model. We then de- 
scribe the simulation methods used and present results 
of simulations of one step of the DNA walker. The mo- 
tor consists of a single DNA strand (cargo) which moves 
from one complementary strand (stator) to the next 
through toehold-mediated strand displacement. We study 
the process for different distances between stators, dif- 
ferent strengths of the attachments of the stators to 
a surface and for different lengths of the toehold, par- 
ticularly focusing on the consequences of inter-stator 
distance. 

2 Burnt-bridges DNA motor 

The burnt-bridges DNA motor studied in this work is a 
system that produces autonomous, unidirectional mo- 
tion of a single-stranded cargo along a track of single- 
stranded stators. The motor was experimentally realised 
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Fig. 1 (a) A schematic illustration of the stepping process 
of a burnt-bridges motor: (i) The cargo (red) is attached to 
the first stator (blue), which is a complementary DNA strand, 
(ii) The nicking enzyme binds to its recognition sequence, and 
catalyses the hydrolysis of the backbone, (iii) The nicking en- 
zyme and stator fragment have dissociated, (iv) The exposed 
toehold of the cargo binds to the second stator (green), (v) 
The second stator fully displaces the first, and the motor com- 
pletes a step, (b) Decision making at a junction. The cargo 
could step to either the yellow or green stator. By altering or 
blocking the toeholds, one direction or the other can be made 
preferable. 
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Fig. 2 Illustration of the structure of nucleotides in the 
coarse-grained model. The figure shows two nucleotides, with 
a spring between backbone sites to indicate the connectivity 
of bases within a strand (maintained by a finitely-extensible 
spring in the model). Each nucleotide possesses base and 
backbone excluded volume sites, a hydrogen bonding site and 
a stacking site. Red spheres indicate backbone and base ex- 
cluded volume sites - the size of the spheres shows the range of 
the excluded volume interaction for each site. Green and blue 
spheres are centred on the stacking and hydrogen-bonding 
sites - the size indicates the distance at which the stacking 
and hydrogen-bonding interactions are maximised (when the 
surface of spheres from different nucleotides come into con- 
tact). Note that the interaction potentials are dependent on 
the relative orientation of nucleotides, as well as the location 
of the interaction sites. 



in Refs. [24l[30l[35| . and the stepping process of the mo- 
tor is schematically illustrated in Figure [TJ The sta- 
tors can be attached to a DNA duplex or a DNA 
origami surface [501155] . 

Initially, the cargo strand is attached to the first 
stator. A nicking enzyme (N.B bvC lb) is present in the 
solution [63] . These enzymes can bind specifically to 
a certain sequence of double stranded DNA present in 
the stator/cargo duplex, and cut the backbone of the 
stator strand a short distance from the 5' end of the 
stator (6 bases in Refs. [301133, and 8 in Ref. [M])- The 
binding of the shorter stator fragment is unstable at ex- 
periment conditions, and it tends to detach, revealing a 
short toehold on the cargo. The next available stator is 
positioned close to the first stator (around 7 nm in Ref. 
[23] and 6 nm in Refs. [501135] ) and the exposed toehold 
can bind to the next stator. Strand displacement can 
then occur, allowing the cargo to replace bonds to the 
first stator with bonds to the second stator. 

Once the displacement process is complete, the cargo 
is totally detached from the first stator and fully bound 
to the second one. The stepping process can now be re- 
peated, with the stator to which the cargo is attached 
being cut again by the enzyme and the cargo making 
a step to the next stator in line. The backward step 
is now highly improbable, as the preceding stator has 
been nicked and has fewer complementary bases with 
the cargo strand. Used stators therefore get disabled 
as the cargo travels along the track, leading to the de- 



scription 'burnt-bridges'. Directional motion is possible 
because the walker's motion catalyses the hydrolysis 
of the stator's backbone, a free-energetically favourable 
process. 

In the original experiment [24] . three stators were 
attached to a double-stranded DNA track. Fluorescence 
was used to demonstrate that the cargo stepped along 
the track, visiting the stators in order. Stators were 
later attached to a DNA origami surface [SOj, and the 
cargo was observed to move along a 17-stator track by 
atomic force microscopy. Motion along an 8-site track 
was also observed via fluorescence, and was strongly 
suppressed by the removal of a stator. Recently, mo- 
tors have been designed that can choose a pathway at 
a junction based on information either carried by the 
motor itself or provided externally [35] . 

The suggested future applications of DNA walkers 
such as the burnt-bridges motor include a programmable 
chemical synthesis and a molecular realization of a Tur- 
ing machine |64) . 



3 Coarse-grained DNA model 

We use the coarse-grained DNA model introduced in 
Refs. [58] and [61]. It represents a DNA strand as a 
string of nucleotides, each of them a rigid body with 
several interaction sites. The structure of the nucleotides 
is illustrated in Figure [2] Inter-nucleotide interactions 
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Fig. 3 Three strands in a double iielical state with schematic 
illustration of the stabilising interactions. Excluded volume 
interactions, which are not shown, prevent the collapse of the 
system into a small volume, although nucleotides rarely ex- 
perience it directly when in the double-helical state. In this 
image, and others in this article, spheres represent the back- 
bone sites and ellipsoids the base sites. Nucleotides in the 
same strand are coloured identically, and the connectivity is 
illustrated by thin cylindrical tubes. This format of illustra- 
tion is found to be clearer than the more schematic represen- 
tation in Fig. [21 

are expressed through the system's potential 

VdNA (^^-b- + ^«tack + K'xc) + 

+ (^HB+Kr.st. +Kxc + K=x.st.), (1) 

where the first sum is taken over all nucleotides that are 
adjacent along the backbone of a strand and the second 
sum is taken over all other pairs of nucleotides. The 
nucleotides that participate in these interactions within 
the duplex state are schematically shown in Figure |31 
The detailed description of the interactions along with 
their functional form is available in Ref. [53] • Below we 
provide a basic description of each term. 

The backbone potential Vb.b. is an isotropic (finitely- 
extensible) spring potential that mimics the covalcnt 
bonds in the backbone by constraining the distance 
between the backbone sites of the neighbouring nu- 
cleotides within a strand. This constraint is represented 
through the cylindrical tubes connecting the nucleotides 
within a certain strand in our figures, such as Figure [31 

Nucleotides also interact through excluded volume 
potentials associated with the backbone and base of 
each nucleotide, V^^c ^^'^ Vc^c- The interaction radius 
of the excluded volume potentials is shown in Figure 
[21 Hydrogen bonding interactions Vhb occur between 
anti-aligned complementary bases, and stacking inter- 
actions V^tack drive the strand to form helical stacks 
of bases. Together, these two interactions are the pri- 



mary driving force for duplex formation. Cross stacking 
Vcr.st. is an additional secondary stabilising interaction 
within a duplex. Coaxial stacking Vcx.st. captures the 
stacking interaction between nucleotides that are not 
adjacent along the backbone of a strand. VhB; K;tack, 
Vcr.st. and Vcx.st. are all anisotropic and depend both on 
the mutual orientations of nucleotides and the distances 
between the interaction sites. 

The model was fitted at salt concentration [Na^] = 
0.5M, where electrostatic properties are strongly screened. 
Hence, for simplicity, repulsive electrostatic interactions 
are simply incorporated into the backbone excluded vol- 
ume. Most DNA nanotechnology is indeed carried out 
at high salt concentration to suppress repulsion between 
strands. 

Although the model has sequence-specific binding 
(i.e., the hydrogen bonding interaction is negative only 
for Watson-Crick base pairs A-T and C-G), it does not 
otherwise distinguish between different base types in 
terms of interaction strengths. The model has recently 
been extended to include sequence-dependent interac- 
tions for base pairing and stacking in Ref. j60j . In this 
work, we use the "average-base" version of the model, 
since generic features of the system are easier to resolve 
without sequence-dependent complications. 

The model captures the structural, mechanical and 
thermodynamical properties of ssDNA and dsDNA with 
average AT and GC content in the sequence [55]. The 
model was fitted to reproduce melting temperatures of 
short duplexes as predicted by SantaLucia's nearest- 
neighbour model |65II36| for an average sequence. The 
single strands in the model are relatively flexible, allow- 
ing the formation of single-stranded hairpins as well as 
double-stranded helices. 

The model's physical representation of DNA prop- 
erties has allowed it to be successfully used for the study 
of a wide range of DNA systems. In terms of nanotech- 
nology, DNA nano tweezers [21] , kissing hairpins [5^150] 
and a two-footed DNA walker [59] have been simulated. 
From a more biophysical perspective, the nematic or- 
dering transition of dense solutions of short duplexes 
[S5] , DNA overstretching [57] and the extrusion of cru- 
ciform structures under torsion |68| have also been mod- 
elled. The physically reasonable behaviour in all cases 
gives us confidence in applying the model to the burnt- 
bridges motor. 

4 Simulating the burnt-bridges motor 

4.1 System 

Our simulated system consists of three DNA strands: 
the first stator, the cargo strand and the second stator. 
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Fig. 4 Typical configurations that are sampled in simula- 
tions of the stepping process, with the distance between the 
stators 7.1nm and attachment spring constant 131pNnm^^. 
The yellow spheres inside red circles indicate the position of 
the points to which the 3' ends of the stators are attached by 
spring potentials. Double-headed arrows indicate transitions 
which are expected to be reversible under the system condi- 
tions, (a) The cargo (red) strand is attached to the first stator 
(blue), with a 6-base toehold exposed after the nicking of the 
first stator. (b) The cargo's toehold makes 6 bonds with the 
second stator (green), (c) Displacement is nearly complete - 
the cargo has two bonds with the first stator and nineteen 
with the second stator. (d) Displacement is complete: the 
cargo is fully bound to the second stator. 



The representation of these strands by the model is il- 
lustrated in Fig. 01 We simulate the process by which 
the motor steps to the next stator, stages (iii)-(v) in 
Fig. [TJ The first stator has six bases fewer than the sec- 
ond stator, emulating the state shown in Fig. [T](a)-(iii) 
after the small fragment of the stator has dissociated 
after nicking. We model the attachment of the stators 
to a surface via a spring potential, acting on the first 3' 
base of the stators: 



Kpring(fc, ri,r2) 



(2) 



where ri and r2 are the positions of the centre of mass 
of the first 3' nucleotide of the first and second stator re- 
spectively and r° and are the positions to which they 
are attached by the springs. This potential energy is in- 
cluded in the total potential energy of the system. We 
define a set of coordinates such that the stator attach- 
ment points lie in the z = plane, and are separated 
by a distance d in the x direction. We will compare 
three different values of d: 3.3 nm, 7.1 nm and 9.4nm. 
7.1 nm approximately corresponds to the distance used 
in Ref. [53] and we chose 3.3 nm and 9.4 nm to test 
shorter and longer distances respectively. 



In most simulations, we use the spring constant k = 
131pNnm~-'^, chosen so that the variance of the dis- 
tance between the attached nucleotides of the first and 
second stators. (|ri — r2p) — (|ri — r2|)^, is approxi- 
mately equal to the variance of the distance between 
two nucleotides that are 1 1 base pairs away on a strand 
in DNA duplex as simulated with our model at tem- 
perature 37 °C. Choosing the spring constant in this 
manner means that its magnitude is physically sensible 
for a DNA-based system; we will, however, consider the 
consequences of varying it. 

To further mimic the presence of the DNA origami 
substrate, we forbid the cargo and stator strands from 
crossing the z = plane. To achieve this we introduce 
an additional potential 



V 



repulsion 



%z,2, if < 







if 2, > 



(3) 



which acts on the positions of centre of mass = 
{xi, yi, Zi) of all nucleotides i in the simulation. We set 
kr to 1142 pNnm^^, which is sufficiently large to effec- 
tively prevent nucleotides from crossing the z = plane. 
The total potential energy of the simulated system is 
then V{k, d) = Vdna + VspHng 

our simulations are done at temperature 37 °C, the tem- 
perature used in experiment. 

We use the following sequences for the first stator, 
second stator and cargo, respectively: 

5 ' -TCAGCCCAACTAACATTTTA-3 ' 

5 ' -GGAACCTCAGCCCAACTAACATTTTA-3 ' 

5 ' -CGATGTTAGTTGGGCTGAGGTTCC-3 ' 

which correspond to the sequences used in Ref. [30] . In 
the experiments, there is an additional 20-base segment 
at the 5' end of the stator that was used to bind a 
blocking strand that was displaced before the beginning 
of the stepping measurements. We do not include this 
segment in our simulation. The last 4 bases at the 3' 
end of the stators are not complementary to the cargo 
strand, acting as a flexible linker with the surface. Note 
that the two bases at the 5' end of the cargo strand 
are not complementary to the stator strands: in the 
experiment [30] this dangling end was used to attach a 
fluorophore to the system for tracking. 



4.2 Simulation methods 

Unbiased Monte Carlo methods sample different eon- 
figurations of the system with a relative probability 
given by the physically appropriate Boltzmann distri- 
bution (exp(— V^/fcaT"), with T being the temperature. 
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/cb the Boltzmann constant and V the potential energy 
of the configuration) . Properties of a system at a certain 
temperature can then be measured by averaging over 
the configurations obtained. We use the Virtual Move 
Monte Carlo (VMMC) algorithm [691170] (specifically, 
the variant in the appendix of Ref. [7D]). VMMC is an 
approach in which new configurations are proposed by 
moving dynamically selected clusters of particles (nu- 
cleotides in our case), accelerating the equilibration of 
strongly-interacting systems relative to simple imple- 
mentations in which only a single particle is moved 
at each step. The proposed moves used in our simu- 
lations are either rotation around a randomly chosen 
nucleotide's backbone site or linear translation. 

In our simulations we compute the equilibrium free 
energy of the system as a function of the number of base 
pairs between the cargo and the first stator {bic) and 
the cargo and the second stator (62c)- We define a base 
pair as being formed if the hydrogen-bonding energy 
(Vhb in Eq. ([ij) between the two bases is more neg- 
ative than — 0.97fcBT, about 15% of typical hydrogen- 
bonding energies in our model. The free energy of a 
state is related to the probability P that the system is 
found in such a state during simulations by 



F (bic, b2c) /kBT = - log P (5i„ 62c) + logPo 



(4) 



where Pq is an arbitrary normalisation. A relatively 
high free energy thus corresponds to a relatively un- 
likely state. Free-energy landscapes do not completely 
determine the kinetic properties of a system, but they 
are indicative of how the system responds to certain 
changes of parameters. For example, the rate at which 
a process occurs is often limited by the need to pass 
through a high free-energy (improbable) state. Raising 
(or lowering) this barrier by perturbing the system gen- 
erally results in an exponential decrease (or increase) in 
the rate of the process. 

Sampling all relevant states can be difficult even 
with an efficient algorithm such as VMMC, as (meta)stable 
states are often separated by high free-energy barriers 
that are difficult to cross. Simulations can therefore get 
stuck for large portions of the simulation time in a lo- 
cal free-energy minimum. To overcome this problem, we 
use the Umbrella Sampling method ^T] . Instead of sam- 
pling from the Boltzmann distribution, we sample from 
configurations using the weight w(&ic, &2c) exp(— IZ/fcer), 
where u'(6ic, 62c) is an arbitrary biasing potential. w{bic, &2( 
can be chosen to raise the probability of visiting unlikely 
transition intermediates, thereby accelerating equilibra- 
tion. To obtain an unbiased estimate of the free energy, 
one must correct for the applied bias: 



i^(&ic,M AbT = - log 



ftiascd {bic, 62c) 
w{bic,b2c) 



logPo, 



(5) 



where Pbiascd (&ic, &2c) is the probability with which states 
appear in the biased simulation. In this work, we run 
several simulations to determine w(&ic,&2c) by hand to 
enhance sampling of the relevant states. Once a sat- 
isfactory w{bic,b2c) are obtained, long simulations are 
run to collect the final data. 

It is convenient in our case to split the Umbrella 
Sampling protocol into two windows. In the first win- 
dow we study attachment of the cargo's toehold to 
the second stator, and in the second we consider the 
displacement process. To do this, we restrict the sys- 
tem to bic > 14 and &2c < 8 in the first case and 
14 < bic + b^c < 23, bic > 1 and b2c > 5 in the sec- 
ond. The windows are then combined using the overlap 
between the two with the weighted histogram analysis 
method [7^ . 

We note that we do not sample all possible val- 
ues of bic and 62c using these two windows. To do so 
would have been computationally costly without be- 
ing likely to provide much insight. In particular, we 
do not consider the breaking of many base pairs of the 
first stator/cargo duplex unless the cargo is attached to 
the second stator. The high free-energy cost of sponta- 
neous breaking of base pairs at this temperature make 
substantial melting extremely unlikely. For the same 
reason, we do not consider displacement intermediates 
with fewer than 14 base pairs involving the cargo. Fi- 
nally, we do not sample the transition to states with 
bic = 0, when the cargo is detached from the first sta- 
tor. To do so would be especially difficult as it would 
require simulating the reattachment of the cargo to the 
first stator in the absence of a toehold. Once detached 
from the first stator, however, the system is in a configu- 
ration very similar to the initial one, the only difference 
being the extra six base pairs available with the second 
stator. The free-energy change associated with adding 
a base pair to an isolated duplex in our model is known 
to be around 2.3 k-oT at 37°C (this can be seen from the 
slope of the free-energy profile in Fig. [5](b) for the ini- 
tial stages of toehold binding), and so we can roughly 
estimate the free energies of the bic — states from 
those of the 62c = states. Most importantly, however, 
only states with both bic > and 62c > are expected 
to be sensitive to the separation of the stators, as only 
in these states is the cargo stretched between both. 

The Umbrella Sampling simulations reported in this 
work involve at least 3 x 10" attempted VMMC moves. 
Half of the moves are rotations (about a random axis 
with an angle drawn from a normal distribution with 
mean zero and standard deviation 0.2 radians), and half 
are translations (through a random vector with a length 
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drawn from a normal distribution with mean zero and 
standard deviation 0.85 A). 



5 Simulation results 

We first present the free-energy landscape, and a pro- 
file along a one-dimensional pathway, for two stators 
7.1 nm apart attached by springs with stiffness k = 
131pN/nm, to illustrate the basic features of the at- 
tachment of the cargo to the second stator and of the 
branch- migration process. We then compare the free- 
energy profiles for different attachment spring constants, 
different distances between the stators, and different 
lengths of the toehold on the first stator. 



5.1 Stators separated by 7.1 nm 

The free-energy landscape as a function of bic and 62c 
is shown in Fig. [5](a), with major features highlighted 
by inspecting the one-dimensional pathway shown in 
Fig.[5](b). The free energy is normalised to be equal to 
for the case when the cargo has no bonds with the 
second stator. The basic features of the landscape are 
the following. 

— There is a rise in free energy associated with the for- 
mation of the initial base pair with the second sta- 
tor, due to the loss of configurational entropy when 
the first contact is formed. 

— The free energy decreases as successive base pairs 
are formed in the toehold, due to the cooperative 
nature of duplex formation (once the first contact is 
formed, successive base pairs are much more likely). 

— As displacement begins (once the seventh base pair 
is formed with the second stator), there is an ini- 
tial rise in free energy, followed by a plateau. This 
initial rise is a generic feature of displacement re- 
sulting from steric interference at the displacement 
interface, and is the subject of a forthcoming paper 

— Later stages of displacement, after around 15 base 
pairs have formed between the cargo and the second 
stator, involve an increase in free energy of around 

4fcBr. 

An unusual feature of the free-energy profile shown 
in Fig. [S](b) is the increase in free energy towards the 
end of displacement. We attribute this to an increase in 
tension within the system. When the first bond between 
the cargo and the second stator is formed, the contact 
point is far away from the nucleotides that are attached 
to the surface (the 3' end of the stators). It is there- 
fore not difficult for the strands to reach each other at 
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Fig. 5 (a) The free-energy landscape of motor stepping as a 
function of bonds between the cargo and the first and second 
stator. The arrow indicates a pathway through the landscape 
that is used to plot (b). (b) The free-energy profile of dis- 
placement, plotted along the one-dimensional path shown in 
(a). A stage on this path is uniquely specified by the number 
of bonds with the second stator. The final point (23 bonds 
with the second stator) is estimated as discussed in the text, 
not measured. 



the contact point. By contrast, when more bonds have 
formed between the cargo and the second stator, the 
contact point is closer to the 3' end of the stators. Even- 
tually, the length of DNA between the contact point 
and the surface attachments gets so short that main- 
taing the structure causes considerable tension, which 
is free-energetically unfavourable and results in the ob- 
served rise in the profile. 

The role of the attachment of stators to a surface can 
be tested by changing the spring constant of the attach- 
ment to the surface. The results of otherwise identical 
simulations with different spring constants are shown 
in Fig. ini Increasing or decreasing the spring constant 
from k = 131 pN/nm by a factor of about 4 has a very 
small effect on the free-energy profile. The reason is 
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Fig. 6 The free-energy profile of motor stepping for various 
strengths of attachment to the substrate, all using a stator 
separation of 7.1nm. The profile is taken along the path illus- 
trated in Fig[5](a). 

that in this range the attachment spring is fairly stiff, 
and the less costly way for the system to come close 
together is to stretch the single-stranded sections. By 
decreasing the spring constant by nearly two orders of 
magnitude, down to /c = 3 pN/nm, we are able to access 
a regime where the attachment springs are sufficiently 
weak that the strands can relax the tension effectively 
by moving the bases at the 3' end of the stators closer 
together rather than stretching the single-stranded sec- 
tions. We stress that the physically relevant regime for 
stators attached to a DNA duplex or to a DNA origami 
is the high-fc one. 

5.2 Varying the distance between stators 

We now compare the stepping of the motor for three 
different distances between the stators: d = 3.3 nm, 
7.1 nm, and 9.4nm. The free-energy profile along the 
one-dimensional pathway indicated in Fig.[5](a) is shown 
in Fig. [7](b) for aU three distances. Fig. [2](a) shows the 
full two-dimensional free-energy landscape for the case 
oi d = 9.4 nm. All the profiles were produced with the 
same spring constant, namely k — 131pNnm^^. The 
free energies have been normalised to zero when there 
are no bonds between the cargo and the second stator. 
One of the typical configurations sampled by our simu- 
lations for distance 9.4 nm between stators is illustrated 
in Fig. El 

The initial part of the free-energy profile, from 
bonds to 6 bonds with the second stator, is nearly iden- 
tical for all three distances considered. As the second 
stator makes more bonds with the cargo strand, we see 
that the increase in free energy is bigger for larger dis- 
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Fig. 7 (a) Two-dimensional free-energy landscape of step- 
ping for d = 9.4 nm. (b) The free-energy profiles of motor 
stepping along the pathway illustrated in Fig[5](a) for three 
different distances d between the attachment points of the 
stators. The illustrated alternative pathway in (a) shows that 
the profile in (b) for d = 9.4 nm probably overstates the dif- 
ficulty of displacement, although it is still far more difficult 
than for d = 7.1 nm. 

tances d between the stators. This effect is consistent 
with our understanding that the rise in free energy is 
associated with increasing tension within the complex 
due to the need to stretch DNA between the surface 
attachment points and the junction. The snapshot in 
Fig. [5] clearly illustrates the tension in the system at 
later stages of displacement for d = 9.4 nm. 

In fact, inspection of Fig. [7](a) shows that the ten- 
sion for the 9.4 nm case is so great that when the cargo 
has only one bond with the first stator, it is thermo- 
dynamically more favourable for the cargo to be bound 
by only 18 or 19 base pairs to the second stator, rather 
than by 21 base pairs as for the other values of d. 
It is therefore highly probable that a typical displace- 
ment pathway would involve the first stator detaching 
when the second stator has significantly fewer than 21 



Simulating a burnt-bridges DNA motor with a coarse-grained DNA model 



9 




Fig. 8 One of the typical configurations sampled in the sim- 
ulation for distance between stators d = 9.4 nm. The cargo 
has three bonds with the first stator and eighteen bonds with 
the second stator. The tension acting on the DNA between 
the attachment points and the displacement interface can be 
clearly seen from the stretched arrangement of the second 
stator (green) close to its attachment point. 




Fig. 9 Illustration of the reduction in constrained bases as 
simulation progresses, (a) Toehold-binding only: the DNA 
along the path X — )> y — i> Z is constrained by the need to 
for the system to from a closed loop X ^ Y ^ Z ^ X (with 
the vector between attachment sites as part of the loop), (b) 
Later on during displacement, fewer bases are constrained 
within the loop. Those that remain, however, are under ten- 
sion due to the need to stretch between attachment sites. 

base pairs with the cargo, along an alternative pathway 
such as that shown in Fig. [3(a). As such, the profile 
along the pathway shown in Fig. [7](b) tends to over- 
state the difficulty in stepping between the two stators 
for d — 9.4 nm, although it is still a much more difficult 
process than for the shorter values of d. 

Interestingly, for the smallest separation of d = 3.3 nm 
there is actually a decrease in free energy with increased 
binding of the cargo to the second stator. We attribute 
this to the fact that the number of bases under con- 
straint actually decreases as displacement proceeds. As 
can be seen from Fig. |9](a), when the second stator is 
bound only to the toehold of the cargo, most of the 
bases in the system are effectively held within a closed 
loop (X — > y — ^ Z — > X in the diagram - this loop 
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Fig. 10 The free-energy profile of motor stepping for two 
different values of cargo toehold length. In both cases, the 
stators are at the same distance of 7.1 nm. For the 6-base 
toehold, the profile is obtained along the path shown in Fig. 
[5](a). For the 9-base toehold, an analogous pathway in which 
the cargo first binds to the second stator by 9 bases and then 
the second stator displaces the first base-by-base is used. Once 
again, the final state is estimated as discussed in the text. 

includes the vector between attachment sites) that re- 
duces their conformational freedom. As displacement 
proceeds, a single-stranded section of the first stator 
that is not constrained by looping is generated, and the 
number of bases in the loop is reduced, as depicted in 
Fig. in](b). All else being equal, transferring bases from 
within the constrained loop to an unconstrained section 
should be a favourable process because of the increase 
in a configurational entropy. However, the extra tension 
felt by the short loop overwhelms this effect for large 
stator separations. 

5.3 Different toehold lengths 

Finally, we compare the free-energy profile obtained at 
d = 7.1 nm and k = 131pNnm~^ with that for an 
identical system, except with a longer toehold (9 bases 
rather than 6) exposed by the nicking of the first stator. 
The free-energy profiles (taken along the same pathway 
as indicated in Fig. [5](a) for the 6-base toehold and 
along an analogous path for the 9-base case) are shown 
in Figure [TUl These profiles have been normalised to 
when the cargo strand is bonded fully by its toehold 
to the second stator, i.e. by 6 and 9 bases respectively. 
As expected, in the case of the 9-base toehold, we ob- 
serve a larger barrier for detachment of the cargo strand 
from the second stator once it is bound by the toehold 
due to the extra base pairs. By contrast, once displace- 
ment is well underway the free-energy profiles do not 
differ by more than 0.3fcBT, which is comparable with 
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the estimated errors. Thus, once displacement has been 
initiated, the free-energy changes are only dependent on 
the distance between the strands and not the toehold 
length. 

5.4 Consequences of free-energy profiles for motor 
operation and track design 

Our data suggest that the tension generated within the 
motor as displacement proceeds has a potentially sig- 
nificant effect on motor operation. In particular, the re- 
sultant rise in free energy will suppress the speed with 
which the second stator can fully displace the first after 
binding to the toehold of the cargo. If this suppression 
is strong enough, the probability that the second sta- 
tor detaches from the toehold rather than completing 
the displacement increases. The slope of the free-energy 
profile changes rapidly with d, being fairly gentle at 
d = 7.1 nm and very significant at d = 9.4 nm. Wc 
would therefore expect its effects on stepping success 
to become noticeable for d ~ 8 or 9nm for a 6-base 
toehold. By contrast, we note that the initial toehold 
contacts have almost equal free energies for the three 
distances studied. So, even though the distances vary, 
this result suggest that initial binding rates are rela- 
tively similar in all cases. 

If stepping to the next stator after the current one 
has been cut is the limiting stage in motor operation, 
then the reduction of this rate would be manifested in 
the overall stepping speed of the motor. However, the 
binding, action and unbinding of the nicking enzyme all 
contribute to the time required to take a step, so (at 
least at low enzyme concentrations and moderate val- 
ues of d) this change may have a negligible effect on the 
overall stepping speed. At sufficiently large d, however, 
increasing d further should have a noticeable effect on 
overall motor speed. This limit was clearly reached in 
Ref. [SO]; when removal of a single stator (estimated 
gap ^ 12 nm) dramatically reduced the overall rate at 
which the cargo reached the end of the track. Our re- 
sults suggest this reduction was primarily due to the 
increased difficulty in completing displacement, rather 
than a reduced rate in making contact between the sec- 
ond stator and the toehold of the cargo. 

Interestingly, the potential reduction in the success 
probability of displacement may be advantageous at 
decision-making junctions such as those illustrated in 
Fig. [11(b). We might encourage the cargo to choose sta- 
tor A by relatively destabilising the toehold of B in 
some manner, such as blocking it with another strand 
|35| or by having a different sequence. Leak currents 
will always exist, however, and even blunt-ended strand 
displacement (with no toehold) can occur [73] . 



The selection ratio of stator A to stator B, (^a/b, 
can be modelled by assuming that the toehold of the 
cargo can initially bind to stator i with a rate 7; , where- 
upon the cargo successfully completes the step to sta- 
tor i with probability fi and detaches with probability 
1 — fi- In this case. 



IaJa 
IbIb' 



(6) 



Assuming A and B are equidistant from the original 
stator, there is limited room to adjust ^a/ib- The fact 
that /i < 1 is important. It means that it will be essen- 
tially impossible to chose between two stators if both 
have sufficiently stable toeholds so that /a and /s ^ 1 
(even if there is a large difference in absolute stability 
between the two). Furthermore, it means ipA/B cannot 
be increased arbitrarily by increasing the stability of 
toehold A^ limiting the maximum signal-to-noise ratio. 

Increasing the failure rate of displacement by ad- 
justing d (for both stators A and B equally) will tend 
to reduce fi, and hence the larger toehold stability will 
be required to achieve /i ~ 1. The junction will then 
be able to distinguish between toeholds that are more 
stable than previously, and will have a higher maximum 
signal-to-noise ratio. We note that if the primary con- 
sequence of increasing d was to reduce the binding rate 
of the toehold to the next stator (rather than the suc- 
cess probability once bound) , it would be impossible to 
improve the efficiency of the junction in this way. 



6 Discussion 

We have briefly introduced a coarse-grained model of 
DNA, and demonstrated that it can be used to investi- 
gate an active DNA nanotechnological device: a unidi- 
rectional molecular motor [^I30[I55] . In particular, we 
have studied the physics of the stepping process from 
one stator to the next, once the first has been cut by 
a nicking enzyme, with emphasis on the dependence of 
this process on the separation of stators. 

We observed that the free-energy profiles of initial 
binding of the cargo's toehold to the second stator are 
fairly insensitive to stator separation. However, as dis- 
placement proceeds, there is a rise in the free energy 
when stators are separated by larger distances, associ- 
ated with the need for ever shorter sections of DNA to 
extend across the gap between attachment points. Such 
a rise will tend to reduce the speed at which the cargo 
completes its step to the next stator once it is bound 
by its toehold. For a large enough rise, there will be 
a reduction in the probability of successful completion 
of a motor's step following initial attachment. These 
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results suggest that the experimentally observed reduc- 
tion in stepping rate when the distance between stators 
is doubled |30j is predominantly due to the increased 
difhculty of completing displacement, rather than a re- 
duced probability for forming an initial contact between 
the next stator and the toehold of the cargo. We argue 
that such a reduction in successful stepping could be 
used to provide higher sensitivity at junctions where 
motors must chose between two adjacent stators. 

We find that the rigidity with which stators are held 
in place can be important. At low stiffnesses, stators 
can move towards each other and relax some of the 
tension generated by displacement. At high stiffnesses, 
the stators' movement is limited and the DNA must 
stretch across the full gap between attachment points 
instead, resulting in the aforementioned rise in free en- 
ergy. These results suggest that the flexibility of the 
surface to which the stators are attached, and the na- 
ture of the attachment, could be significant in deter- 
mining the properties of motor stepping. We estimate, 
however, that anchoring stators such as those studied 
here to a single duplex is enough to qualify as a stiff 
attachment. 

Our results have been obtained with a coarse-grained 
model, and the free-energy landscapes and profiles we 
have measured do not completely determine kinetics, 
which also depends on non-equilibrium effects. Nonethe- 
less, the increase in tension within the motor as dis- 
placement proceeds is based on very general mechani- 
cal and structural properties of DNA that are known 
to be well reproduced by the model. So it is reason- 
able to assume that large stator separation does indeed 
lead to a rise in free energy with displacement progress. 
Furthermore, such a rise makes later displacement in- 
termediates harder to reach and will eventually reduce 
the probability of successful step completion once the 
second stator is bound to the cargo's toehold. For the 
system considered in this work, a separation of around 8 
or 9 nm should be sufficient to demonstrate the effects 
of stator separation. The most important assumption 
that we have made is that the small fragment of the 
first stator and the nicking enzyme tend to dissociate 
before attachment to the next stator can occur. Even 
if this is not the case, however, their infiuence will be 
strongest during toehold attachment and displacement 
initiation, and should not prevent the rise in free energy 
associated with increased tension at the later stages of 
displacement. 

Future work could include explicit simulations of the 
dynamics of motor stepping, although such simulations 
will probably be computationally demanding even for 
a model as simple as this one. It would also be worth- 
while to study the consequences of different sequences 



using the sequence-dependent model, and motifs such 
as internal mismatches, on the operation of the motor. 

The software tools implementing our coarse-grained 
DNA model that were used for producing data for this 



work are freely available at http: / /dna. physics. ox. ac.uk 
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